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Entanglement is one of tlie key resources required for quantum computation[l], so experimentally 

creating and measuring entangled states is of crucial importance in the various physical implementa- 
tions of a quantum computer [2]. In superconducting qubits [3], two-qubit entangled states have been 
demonstrated and used to show violations of Bell's Inequality [4] and to implement simple quantum 
algorithms [5]. Unlike the two-qubit case, however, where all maximally-entangled two-qubit states 
are equivalent up to local changes of basis, three qubits can be entangled in two fundamentally differ- 
ent ways [6], typified by the states |GHZ) = (|000)-h|lll))/\/2 and |W) = (|001)-K|010)-h|100))/v/3. 
Here we demonstrate the operation of three coupled superconducting phase qubits [7] and use them 
to create and measure |GHZ) and |W) states^. The states are fully characterized using quantum 
state tomography [8] and are shown to satisfy entanglement witnesses [9], confirming that they are 
indeed examples of three-qubit entanglement and are not separable into mixtures of two-qubit en- 
tanglement. 



In order to create arbitrary entangled states or perform 
arbitrary computations, a quantum computer must im- 
plement a set of universal gates [1], typically taken to be a 
two-qubit gate such as controlled-NOT (CNOT) plus sin- 
gle qubit rotations[ll]. Alternately, universality is possi- 
ble using a three-qubit gate such as the TofiPoU gate[12- 
14]. Three-qubit gates are also important in such ap- 
plications as quantum error-correction [15] and they can 
simplify some quantum circuits[13]. Because supercon- 
ducting phase qubits can be coupled simply by connect- 
ing them with a capacitor[7], we can design multi-qubit 
interactions that directly generate multi-qubit gates[16], 
rather than building them up from more elementary two- 
qubit gates. For creating the two types of three-qubit en- 
tanglement we employ both approaches, using two-qubit 
gates for the |GHZ) protocol, but using a more efBcient 
entangling protocol for ]W) based on a single three-qubit 
gate. 

The I GHZ) protocol[17, 18] is shown as a quantum cir- 
cuit diagram in Figure la. Starting in the ground state 
1 000), a rotation is applied to qubit A to create the su- 
perposition (jOOO) + |100))/v^. Next, a CNOT gate is 
applied to flip qubit B conditioned on qubit A, resulting 
in the state (jOOO) + |110))/\/2. Finally a second CNOT 
is applied to flip qubit C conditioned on B, resulting in 
the desired state |GHZ). As is typical with quantum cir- 
cuits, this is written in terms of CNOT gates which take 
a simple form in the qubit basis. In our system however, 
a more natural universal gate is the so-called iSWAP 
gate[19] which is generated directly by applying the avail- 
able coupling interaction iJ^f = {hg/2){a^a^ + cFyCfy) 



^Independently, entanglement created between three supercon- 
ducting transmon qubits is being reported in a simultaneous 
publication[10]. 



for time tiswAP = 7r/2g, where g is the coupling strength. 
The I GHZ) protocol can be "recompiled" in terms of this 
gate to obtain the circuit shown in Figure lb. 

The protocol to generate a )W)-state, shown in Figure 
Ic, is based on two features of the )W)-state: it is sym- 
metric with respect to permutations of the qubits, and 
it is a superposition of three states each with one qubit 
excited. Thus, generating the state requires "sharing" a 
single excitation symmetrically among three qubits. This 
is done by first applying a 7r-pulse to qubit B to excite 
it with one photon and create the state jOlO). Then the 
qubits are entangled by turning on an equal interaction 
between all pairs i?int = H-^^ + H^^ + ffj^f for time 
tw = (4/9)tiswAP- The interaction causes the excita- 
tion to be distributed among the qubits, and at time 
the system is left in an equal superposition state, as de- 
sired. A final Z-rotation can then be applied to correct 
the phase of qubit B, though this does not affect the en- 
tanglement of the state. This protocol requires only a 
single entangling operation, and the interaction is only 
applied for a short time, shorter even than the character- 
istic time for two-qubit gates in the system. This yields 
a highly efficient state-generation protocol based on the 
multi-qubit gate generated by iJint- 

To allow for future expansion beyond the present work, 
the sample was designed with four qubits, so that the 
coupling network for the desired symmetric coupling be- 
tween all pairs of qubits is as shown in Figure Id (left). 
The design can be simplified by transforming the cou- 
pling network into an equivalent circiiit (right) with each 
qubit coupled capacitively to a central "island". This 
simplified design is easier to lay out symmetrically on 
chip and requires only N capacitors to couple N qubits, 
rather than the N{N — l)/2 capacitors in the complete 
network. 

Figure 2a shows the complete schematic of the device 
with four phase qubits connected by the capacitive island 
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FIG. 1: Protocols for generating entangled states, a, Quan- 
tum circuit for generating |GHZ) using CNOT gates, b, 
Quantum circuit for |GHZ) that has been "recompiled" to 
use iSWAP gates, which are directly generated by capacitive 
coupling in the phase qubit. These two circuits are not fully 
equivalent, but they both produce |GHZ) when operating on 
the ground state as input, c, Circuit to generate |W) using 
a single entangling step with simultaneous coupling between 
all three qubits. The entangling operation is turned on for a 
time ~ (4/9)fiswAP where tiswAP is the time needed to 
complete an iSWAP gate between two qubits. d, Capacitive 
coupling network to achieve symmetric coupling between all 
pairs of qubits (left), and simplified equivalent circuit using 
coupling to a central island (right). The complete network 
on the left requires six capacitors, and the coupling strength 
g is proportional to the qubit-qubit capacitance Ca- In the 
equivalent circuit on the right, the same coupling strength is 
attained by scaling the capacitors to Cc = 4Ca , but now only 
four capacitors are required and the circuit can be easily laid 
out symmetrically on a chip. 



coupler. Each qubit is individually controlled by a bias 
coil which sets the operating flux bias and also carries mi- 
crowave pulses for manipulating and measuring the qubit 
state. In addition, each qubit is coupled to an on-chip su- 
perconducting quantum interference device (SQUID) for 
state readout. Figure 2b shows a micrograph of the fab- 
ricated device, made from aluminum films on sapphire 
substrate with Al/A10a;/Al Josephson junctions. The 
completed device is mounted in a superconducting alu- 
minum sample holder and cooled in a dilution refrigerator 
to ~ 25 mK. Bringup and calibration of the multiqubit 
device are similar to previous works[20, 21]. Although 
the coupling capacitors are fixed, the effective interaction 
can be controlled by tuning the qubits into resonance at 
/b = 6.55 GHz (coupling "on") or by detuning A and 
C to ±250 MHz (coupHng "off") [22]. The measured cou- 
pling strengths were found to be within 5% of 12.5 MHz 
for each pair of qubits. Importantly, all qubits can be 
brought into resonance simultaneously, as required for 
the |W) protocol, or two qubits can be tuned into reso- 
nance with the third detuned, as required for the iSWAP 
gates in the |GHZ) protocol. 

Capacitive coupling as used here is simple and well- 
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FIG. 2: Device description and operation, a. Schematic of 
coupled-qubit circuit. Each qubit is controlled individually 
by a flux bias line which sets the DC operating point, and 
provides quasi-DC pulses for tuning the qubits in and out of 
resonance and AC (microwave) control signals for qubit rota- 
tions. In addition, each qubit is coupled to a superconducting 
quantum interference device (SQUID) for readout of the qubit 
state. The qubits are capacitively coupled to the central is- 
land, which results in symmetric coupling between all pairs 
of qubits. b, Photomicrograph of the sample, fabricated with 
aluminum (light areas) on sapphire substrate (dark areas). 
The coupler is the cross-shaped structure in the center, and 
the simplicity of this design is evident in the straightforward 
correspondence between the schematic and the completed de- 
vice. The entire sample is mounted in a superconducting alu- 
minum box and cooled to 25 mK in a dilution refrigerator. 



understood but is subject to measurement crosstalk [4, 7], 
which can cause, for example, a state |001) to be er- 
roneously read out as |011), |101) or even |111). This 
crosstalk affects measured probabilities of all excited- 
state populations, however it has no effect on the "null- 
result" probability of measuring |000), since crosstalk can 
only act if at least one qubit is excited. By measuring 
various subsets of qubits and recording the null-result 
probability for each subset, we are able to reconstruct 
the combined state occupation probabilities without any 
effect from measurement crosstalk^. 

Figure 3 shows the time-evolution of the state occupa- 



See supplementary information. 
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tion probabilities during the entangling protocols, using 
crosstalk-free measurement. In the |W) protocol (Figure 
3a) , one qubit is excited and then the symmetric interac- 
tion between all pairs of qubits is used to distribute that 
excitation among all three, as described above. When the 
interaction time is chosen properly, the system reaches an 
equal superposition, and subsequently stays there while 
the interaction is off (Figure 3b). Because this one ex- 
citation is swapped among the various qubits, the state 
evolution during this protocol is clearly visible in the oc- 
cupation probabilities as they evolve in time. 

Figure 3c shows the state occupation probabilities dur- 
ing the I GHZ) protocol, plotted in segments correspond- 
ing to the stages of the protocol as indicated. The ini- 
tial rotations create an equal superposition of all qubit 
states, with all probabilities converging on 1/8. The ef- 
fect of the two iSWAP gates is then primarily to adjust 
the phases of the various components of the superposi- 
tions, so that in the final rotation constructive interfer- 
ence causes |000) and |111) to be populated, while all 
other states are depopulated. The occupation probabili- 
ties behave as expected, but most of the state evolution 
is hidden in the phase information not captured by these 
probability measurements. 

To fully characterize the quantum states created by the 
entangling protocols, including the phase information, 
we perform Quantum State Tomography (QST) by ap- 
plying various combinations of single-qubit rotations be- 
fore measurement. The density matrix is extracted from 
the measured data using maximum likelihood estimation 
(MLE) to find the state that best fits the data while also 
satisfying the physicality constraints that it be Hermi- 
tian positive semi-definite with unit trace. Using this 
procedure, we extract pw and pghz, shown respectively 
in Figures 4a and 4b. Comparing the measured states 
with theory, we find fidehties Fw = (W| pw |W) = 0.78 
and Fghz = (GHZ| pcHZ |GHZ) = 0.62. 

To understand the significance of the measured fideli- 
ties, we compare these results to entanglement witness 
operators that detect three-qubit entanglement. Thrcc- 
qubit entanglement is witnessed[9] for the |W)-state pro- 
vided that Fw > 2/3, and for the |GHZ)-state pro- 
vided that Fqhz > 1/2. Both inequalities are satis- 
fied by the respective measured density matrices, indi- 
cating that they are genuine three-qubit entangled states 
that cannot be decomposed into mixtures of separable 
states. In addition, pcHZ is found to violate the Mermin- 
Bell inequality [23], as predicted by quantum mechanics 
but disallowed by the classical assumptions of local real- 
ity^. The violation is not loophole- free due to use of the 
crosstalk-free measurement protocol rather than a simul- 
taneous measurement protocol [4], but it is nonetheless 
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FIG. 3: Generation of entangled states in the time domain. 
In each row, the left panel shows the pulse sequence with 
time on the horizontal axis and qubit frequency on the verti- 
cal axis. The right panel shows the measured state occupation 
probabilities as a function of time during each sequence, a. 
To characterize the three-qubit interaction, all qubits are ini- 
tially detuned and qubit B is excited with a 7r-pulse. The 
qubits are then tuned into resonance to turn on the interac- 
tion for some time, then detuned and measured. During the 
interaction, the excitation from qubit B (|010)) is swapped to 
qubits A and C (|100) and |001)), then back again. Proba- 
bilities Pioo and Pooi are nearly equal throughout the entire 
sequence, indicating that the coupling is nearly symmetric, as 
desired. At the first crossing point where the three probabil- 
ities are equal, the system is in a |W)-state. b, The coupling 
is turned on until the crossing point is reached and then the 
qubits are detuned, leaving the system in a |W)-like state, 
up to phase rotations due to the detunings. The small resid- 
ual oscillations visible after the qubits have been detuned are 
due to interactions with |2) and higher energy levels of the 
phase qubits. c. The jGHZ) sequence is a direct translation 
of the circuit shown in Figure lb, with the iSWAP gates im- 
plemented by tuning the qubits pairwise into resonance for 
time fiswAP ~ 40 ns. The occupation probabilities on the 
right are plotted versus time in each marked stage of the se- 
quence. After creating the initial superposition (1), the two 
iSWAP gates change the phases of the various components 
of the state, with little effect on the populations (1-2, 2-3). 
During the final rotation, constructive interference populates 
|000) and | 111), while destructive interference depopulates the 
other states. For an ideal |GHZ) state, the probabilities Pooo 
and Pill should approach 50%, though in the experiment this 
level is reduced due to the effects of decoherence and errors 
discussed in the text. 



an indicator of genuine three-qubit entanglement. 

The lower fidelity of jGHZ) compared to |W) is due 
to two main factors: first, the jGHZ) sequence is longer 
because of the two iSWAP gates; the sequence length 
is a substantial fraction of the dephasing time T2 of 
the qubits, which is particularly harmful because the 
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FIG. 4: Quantum state tomography of |GHZ) and |W). At left, the real parts of the measured density matrices pw (a) and 
Pghz (b) are shown in the bar plots. For both states the theoretical density matrix has vanishing imaginary part, and the 
measured imaginary parts (not shown) are also found to be small, with |Impw| < 0.03 and |Im/9GHz| < 0.10 respectively. At 
right, the Pauli set or generalized Stokes parameters are plotted for pw (c) and pghz (d). The bars show expectation values of 
combinations of Pauli operators on one, two and three qubits, with theory in gray and experiment overlayed in color. The same 
state information is contained in both representations, but the Pauli sets clearly show the differences between |W)-type and 
|GHZ)-type entanglement. In addition to the three-qubit correlation terms, the |W)-state has two-qubit correlations because 
tracing out one qubit from a |W)-state still leaves the others partially entangled. The fidelity is Fw = 0.78. For |GHZ), the 
two-qubit correlations other than the trivial ZZ-type are absent because tracing out one qubit leaves the others in a completely 
mixed state. The fidelity is -Fqhz = 0.62 and the state is also found to violate the Mermin-Bell inequality [23, 24]. 



sequence relies on precise phase adjustment and inter- 
ference to populate 1 000) and |111) while depopulating 
all other states. Longer coherence times would improve 
this, as would stronger coupling to reduce the gate time. 
Second, the presence of |2) and higher levels and the rel- 
atively small nonlinearity of the phase qubit cause errors 
due to transitions into higher excited states, for example 
1 110) — > |200). These transitions can be ignored in the 
|W) protocol since they are inaccessible with only one ex- 
citation in the system, but they cause errors in the |GHZ) 
protocol since all qubit states are populated, including 
those with multiple excitations. The effect of higher lev- 
els becomes particularly complicated in this experiment 
when using fixed capacitive coupling with detuning to 
turn off the interaction, due to spectral crowding from 
the higher qubit levels. This highlights the need to re- 
place frequency detuning with tunable coupling schemes, 
which are currently an active area of research. 

In conclusion, we used superconducting phase qubits 
to generate both types of three-qubit entangled states, 
namely |GHZ) and |W). In both cases, the created states 
violate entanglement witnesses that rule out biseparabil- 
ity, showing that these are genuine three-qubit entangled 
states. This ability to couple three qubits and create 
entangled states with qualitatively different types of en- 
tanglement represents an important step toward scalable 



quantum information processing with superconducting 
devices. 
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MATERIALS AND METHODS 

The phase qubits in this device were designed to have 
critical current /q ~ 2 /zA, capacitance C ~ 1 pF, and 
inductance L w 720 pH. The coupling capacitance was 
Cc ~ 15 pF, chosen to give coupling strength 2g/2TT « 
15 MHz at a qubit frequency of 6 GHz. The fabrica- 
tion process is the same as has been used in previous 
experiments[l], using a sapphire substrate with super- 
conducting Al films, AlOx tunnel junction, and a-Si:H 
dielectric for the qubit and SQUID shunt capacitors and 
for wiring crossovers. Each qubit is controlled with a sin- 
gle line which provides DC flux bias, quasi-DC detuning 
and measurement pulses, and microwave state-rotation 
pulses. Each qubit is coupled to a three-junction mea- 
surement SQUID for readout. 

The microwave control signals arc produced by a cus- 
tom microwave arbitrary waveform generator (AWG), 
which has been described in the supplementary infor- 
mation to previous work[2]. For single-qubit rotations 
we use 6 ns full-width at half-max (FHWM) gaussian 
pulses. To reduce errors in these operations due to the 
presence of the 1 2) -state and the small nonlinearity of the 
phase qubit, the pulses include a quadrature modulation 
correction inspired by the scheme known as Derivative 
Removal by Adiabatic Gate (DRAG) [3]. 



BRINGUP AND CALIBRATION 

To operate the sample, we begin by characterizing the 
operation of the individual qubits. Prom basic single- 
qubit experiments, we find relaxation times and spin-echo 
dephasing times for the qubits, as shown in Table I. The 
single-qubit characterizations were performed with the 
qubits detuned by ±250 MHz as indicated in the table, 
so that the coupling interaction is off since the detuning 
A is much larger than the coupling strength A/gi w 20. 

Turning the coupling interaction on requires tuning the 



qubit 


Ti (ns) 


Techo (ns) 


/lo (GHz) 


A 


460 


270 


6.2995 


B 


460 


300 


6.5506 


C 


450 


390 


6.7988 



TABLE I: Single-qubit parameters. 



qubits 


iiSWAP (ns) 


2.g/27r (MHz) 


AB 


40.3 


12.4 


AC 


40.9 


12.2 


BC 


38.8 


12.9 



TABLE II: Qubit-qubit coupling parameters. 

qubits into resonance with each other for some time. The 
required detuning pulses are calibrated pairwise between 
the qubits; for example, to calibrate the detuning pulse 
to couple A with B, we first excite qubit B with a 7r-pulse 
and then adjust the amplitude and length of the detuning 
pulse on A to maximize the transfer of this excitation to 
A. This pulse creates an iSWAP gate, as needed for the 
I GHZ) protocol. The process is repeated for each pair 
of qubits, giving swap times and coupling strengths as 
shown in Table II. The coupling strengths are within 5% 
of each other and also quite close to the design value of 
15 MHz. For the |W)-state protocol the interaction time 
is % = (4/9)*iSWAP « 18 ns. 

CROSSTALK-FREE MEASUREMENT 
PROTOCOL 

The capacitive coupling used in this sample is subject 
to measurement crosstalk, in which the measurement of 
one qubit as being in state |1) can radiate energy into the 
circuit and cause other qubits to switch. With only two 
coupled qubits, crosstalk can cause errors |01) — >■ |11) and 
1 10) — > 1 11). In previous work with coupled qubits [4], we 
have measured the probabilities of these two errors and 
then simply corrected the qubit measurement results to 
account for this crosstalk. However with throe qubits 
all coupled together, the effect of crosstalk is much more 
complicated because there are more possible errors, for 
example |001) |011), |001) ^ |101), |001) |111) and 
so on. While all these error probabilities could in princi- 
ple be measured and corrected for, we instead choose to 
measure in a way that is entirely insensitive to measure- 
ment crosstalk, thus eliminating the need for any correc- 
tion or even any characterization of the crosstalk errors. 

The essential element of this scheme is the fact that 
crosstalk does not affect the state |000) since no qubits 
are excited. Hence when all the qubits are mea- 
sured to find the probabilities PqoOj ^ooi, •••An of var- 
ious outcomes, only the "null-result" probability Pqoo 



2 



of no qubits excited requires no crosstalk correction. 
Now, consider the case when qubits A and B are mea- 
sured but qubit C is not. We obtain four probabil- 
ities Pqox, Poix, Pwx, Piix: where the subscript x indi- 
cates that in each case we have no information about 
the unmeasured qubit C. As before, only the null-result 
probability Ppoa: is unaffected by crosstalk, but note that 
-Pooa; = -Pooo + -Pooi because there are two possibilities for 
the state of the unmeasured qubit C. The two probabili- 
ties Poox and Pooo can be measured without crosstalk as 
just described; from these the third probability Pqoi can 
be determined. 

Continuing in this manner, we can reconstruct 
the complete set of occupation probabilities without 
crosstalk by repeating the experiment 2^ — 1 = 7 
times, each time measuring only a certain subset of 
the qubits and recording only the null-result probabil- 
ity for that subset (for the degenerate case in which 
no qubits are measured, we have Pxxx = !)• This 
gives the following set of measured null-results PnuU = 
{Pqoo, Poox, PqxOt Pqxx, Pxoo, PxOxi PxxO, Pxxx)'^ 1 which is 
related to the desired set of occupation probabilities 

P = (PooO) Pool) PoiO; Poll) Pioo, Pioi, PiiO) Piii)^ ac- 
cording to 

/lOOOOOOOX 

1 1 
1 1 
11110000 
Pnull - 10001000 

110 110 
10 10 10 10 

\ 1 1 1 1 1 1 1 1 / 

By inverting this equation, we thus obtain the occupation 
probabilities in a way that is completely insensitive to 

measurement crosstalk. 

Even with one qubit and no crosstalk, the measure- 
ment of the qubit state has finite fidelity, due to stray 
tunneling of the |0)-state and relaxation of the |l)-state 
during the measurement pulse. Measured results can be 
corrected to account for this finite fidelity, and this cor- 
rection is simpler than the crosstalk correction as the 
fidelity errors act independently on each qubit. Hence 
the inulti-qubit correction is just like single qubit case. 
For each qubit, we determine the fidelity of |0)-state mea- 
surement /o and the error probability of |l)-state mea- 
surement ei = 1 — /i from the so-called measurement 
"s-curves" [5] . Taking these into account. Equation SI 
relating the measured null results to the actual occupa- 
tion probabilities becomes 

(S2) 



Inverting this equation gives the occupation probabili- 
ties with single-qubit measurement fidelity taken into ac- 
count. 

The final subtlety in this measurement process is that 
because each element of PnuU is measured in separate 
repetitions of the experiment, each is subject to indepen- 
dent statistical noise. As a result, when Equation S2 is 
inverted the elements of P may not be nonnegative and 
they may not sum to unity, as required for the set of 
probabilities. We thus use maximum likelihood estima- 
tion (MLE)[6] to enforce these constraints and find the 
probabilities P which would give the measured results 
PnuU with the highest probability. 

QUANTUM STATE TOMOGRAPHY 

To reconstruct the density matrix of the quantum state 
after the entangling protocol, we use a standard imple- 
mentation of Quantum State Tomography [6]. After the 
state is prepared, we apply combinations of the identity 
operation or tt/ 2- rot at ions about X or F to each qubit 
before measuring the probabilities as outlined in the pre- 
vious section. MLE is used to extract a density matrix 
that gives the measured data with highest probability 
while also obeying the constraints on a physical density 
matrix: it must be Hermitian, positive semidefinite (no 
negative eigenvalues) and have unit trace. If we forego 
the MLE procedure and instead use a straightforward 
least-square matrix inversion to perform state tomogra- 
phy, the resulting density matrix is typically Hermitian 
with unit trace, but may have one or two small negative 
eigenvalues (with magnitude on the order of 0.05 times 
the largest eigenvalue), indicating that these are likely 
due to independent statistical noise and finite fidelity of 
the individual measurement results, and not large sys- 
tematic effects. The theoretical and experimental density 
matrices are shown in Figure SI. 

CHARACTERIZING ENTANGLED STATES 
Entanglement Witnesses 

An entanglement witness is an operator that detects 

entanglement in quantum states, either pure or mixed. 
We give here a few basic facts about witnesses, following 
closely the discussion in [7]. An operator W is an entan- 
glement witness if TiiWp) > for all separable states 
p, but Tr{Wp) < for some entangled states p in which 
case we say that the entanglement is "detected" by W. 
Note that such an entanglement witness defines a suffi- 
cient but not necessary condition for p to be entangled, 
and in general there are many possible entanglement wit- 
nesses that each detect some subset of entangled states. 
Some entanglement witnesses are constructed in such a 
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FIG. SI: Quantum state tomography of jW) and jGHZ). a, The real part of pw = |W) (Wj. All nonzero elements are 
equal to 1/3, and all imaginary parts (not shown) are identically zero. The experimental pw real part (b) and imaginary 
part (c) compare nicely with the theoretical prediction, giving fidelity Fw = 0.78 with |Impw| < 0.03. d. The real part 
of Pghz = |GHZ) {GHZ|. All nonzero elements are equal to 1/2, and all imaginary parts (not shown) are identically zero. 
The experimental pghz real part (e) and imaginary part (f) again agree nicely with theory, giving fidelity Fghz ~ 0.62 with 
|ImpGHz| < 0.10. 



way that they can be efficiently measured [8], allowing en- 
tanglement to be detected without performing full quan- 
tum state tomography which requires a number of mea- 
surements that grows exponentially with the number of 
qubits. In the present experiment, however, we perform 
full state tomography on the three-qubit states, as the re- 
quired number of measurements is still quite manageable, 
and then the measured density matrix p can be checked 
against any desired entanglement witness. 

The entanglement witnesses we use all have the form 
W = el — IV') {ip\ where 1-0) is some pure entangled state, 
1 is the identity operator, and e is a constant that de- 
pends on the state \ip). For such a witness, the detection 
condition reads Tr(yVp) = e — {ip\ p \ip) < which then 
becomes = (i/'|p|V') > £ where is the fidehty of 
p with respect to the state \4>)- Hence, entanglement is 
detected if the fidelity is sufficiently high. For the states 
of interest the relevant witnesses are [7] 

Ww = ^1-|W)(W| (S3) 

WghZi = ^1-|GHZ)(GHZ| (S4) 

Wghz, = ^1-|GHZ)(GHZ|. (S5) 

We have changed the names of these witnesses from those 
given in the reference because our emphasis is slightly 
different. For our purposes, Ww is a witness that distin- 
guishes |W)-like states from separable states and likewise 
VVghZi is a witness that distinguishes |GHZ)-like states 



from separable states. Recalling from above that these 
witnesses are satisfied for fidelities greater than e, we 
see that the experimental fidelities Fw = (W| pw |W) = 
0.78 > 2/3 and Fghz = (GHZ| pghz |GHZ) = 0.62 > 1/2 
both satisfy the respective entanglement witness. 

The final witness operator >Vghz2 is a stronger condi- 
tion that distinguishes the class of |GHZ)-like states from 
the class of |W)-like states. The experimental fidelity 
Fghz < 3/4 does not satisfy this witness, so that at least 
according to this criterion it is not possible to separate 
our PGHZ from a convex combination of |W)-like states. 
Recall however that entanglement witnesses always give 
sufficient but not necessary conditions for identifying en- 
tanglement. A better measure of the tripartite entangle- 
ment of Pghz is provided by the Mermin-Bell inequality. 

Mermin-Bell inequality 

The I GHZ) state was first discussed [9] for its very 
strong quantum correlations which could rule out hidden 
variable models with a single measurement, rather than 
by taking many measurements and looking at correla- 
tions among them, as with Bell's Inequality. However, 
this single-measurement violation would require an ideal 
pure state and perfect measurement fidelity, neither of 
which are experimentally feasible. Mermin[10] showed 
inequalities for the realistic mixed-state case that are 
obeyed by hidden variable models, but violated by |GHZ) 
and its generalizations to higher numbers of qubits. He 
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also showed that the potential violation grows exponen- 
tially with the number of qubits. For the three-qubit case 
considered here, classical hidden variable models must 
obey G = {XXX) - {XYY) ~ {YXY) ~ (YYX) < 2, 
while the pure |GHZ) state satisfies G|ghz) = 4. Exper- 
imentally we find GpGHz = 2.076, which violates the in- 
equality and hence rules out a hidden- variable model for 
the measurement correlations of the created state /oghz- 
This violation unambiguously separates the experimental 
state pghz from the class of |W)-like states. 
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